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1. Introduction 

We study entanglement in a one-dimensional <j-deformed valence-bond-solid (VBS) 
state with a spin-S" at each site. This state, which we denote as VBSq(S'), is invariant 
under the action of the generators of the SU g (2) quantum algebra [1,2]. Historically, 
this algebra featured prominently in the analysis of the quantum sine-Gordon model [3] 
and the anisotropic XXZ Hciscnberg chain [4] . Here we use this formalism to introduce 
anisotropy into a spin-S 1 VBS state by deforming the usual underlying SU(2) symmetry 
In a previous analysis of the S = 1 ¥685(5*) state wc used transfer matrix techniques 
to calculate the entanglement spectrum and entropy as functions of the parameter q [5] . 
This investigation and the current one are motivated by two fundamental problems. 
The first is to quantify the effects of anisotropy on the entanglement in VBS states 
[5-9]. Second, we aim to develop effective thermal models [10] that have the same 
spectrum as the density matrix of a subsystem of a VBS state. This last objective is 
motivated by our discovery [5] that the deformation parameter q may be interpreted 
as the temperature of effective models describing the interaction of degrees of freedom 
on the boundary of the subsystem. 
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The undcformcd VBS state is the ground state of the SU(2) symmetric one- 
dimensional AfHeck-Kcnnedy-Lieb-Tasaki (AKLT) model [11,12]. This model has 
nearest-neighbour interactions between integer spin-«S"s. It is described by a 
Hamiltonian of the form H = Y^i hi,i+i- The local Hamiltonian hi^+i is a projector 
onto the subspace spanned by the (5 + 1), (5 + 2), . . ., and 25-multiplcts formed 
by spins at sites i and i + 1 [11-13]. Exact results for the block entanglement 
entropy for the 5 = 1 VBS state were obtained in [14, 15] and for arbitrary integer 
5 in [16,17]. Generalized VBS states with other symmetries, such as SU(iV) [18- 
20], SO(7V) [21,22], Sp(AT) [23], and supersymmetry (SUSY) [24,25] have also been 
considered. Entanglement in the SU(JV) [26,27] and SO (AT) [28] symmetric VBS state 
has been studied, but not in the Sp(iV) and SUSY cases. Since these VBS states 
have matrix product state (MPS) representations [29-33], one can obtain the block 
entanglement in these systems by general transfer matrix techniques [5]. 

The anisotropic g-deformed generalization of the spin-1 AKLT chain was first 
considered in [30,31,34]. The MPS representation of the ground state of the model 
was constructed in [30]. This ground state is separated from excited states by a gap 
[31]. Hence, the spin-spin correlation functions decay exponentially [30, 31]. The exact 
entanglement spectrum of blocks of arbitrary length in the spin-1 VBS g (5) state was 
calculated in [5]. The higher integer spin generalization of the g-deformed AKLT 
model was first proposed in [35,36], where the spin-spin correlation functions were 
calculated. To our knowledge, the entanglement spectrum and entropies of VBS g (5) 
states for arbitrary integer 5 have not yet been evaluated. 

In this paper, we calculate the entanglement spectrum and entropies of g-deformed 
VBS states with arbitrary integer 5. Our analytical approach involves transfer matrix 
methods and the use of (/-deformed Clebsch-Gordan coefficients and 6j symbols. We 
begin by constructing the VBS g (5) state by requiring it to be a ground state of a q- 
deformed spin-5 AKLT Hamiltonian. We shall denote this ground state by the state 
vector |VBS 9 (5)). This state is then partitioned into a block of £ sequential spins 
and the environment E. The density matrix of the whole ground state is therefore 
p = |VBS g (5))(VBS 9 (5)|/(VBS 9 (5)|VBS 9 (5)). We then compute the reduced density 
matrix pi of the block by taking the partial trace of p over the environment, pi = tr# p. 
The details of these derivations are provided in section 2. 

In the double scaling limit of an infinitely long block in an infinitely long chain, 
we are able to exactly diagonalize the reduced density matrix (section 3.1). We then 
use the eigenvalues of this matrix to construct the entanglement spectrum of the block 
(section 3.2). This entanglement spectrum (introduced in [10]) enables us to construct 
an effective Hamiltonian that completely describes the reduced density matrix. The 
eigenvalues of the reduced density matrix are further used to calculate the Renyi and 
von Neumann entanglement entropies 

5 R (a) ee a > 0, (1) 

1 — a 

5 v n = -tr(^lnp^) = lim 5 R (a). (2) 

a— >l 

We thus provide exact measures of entanglement [37-40] in the |VBS g (5)) ground 
state as functions of the deformation parameter q and spin 5. 

We further consider the case of blocks of finite length. We obtain the exact 
eigenvalues of the reduced density matrix in the isotropic case q — 1. With this result, 
we calculate the leading order finite-size corrections to the entanglement spectrum 
and entropies (section 4.1). For a general Q-dcformation, we estimate the eigenvalues 
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of the reduced density matrix in the limit of long but finite blocks by perturbation 
theory (section 4.2). Furthermore, we numerically investigate the properties of the 
reduced density matrix of the spin-2 VBS^S) state. This result allows us to make 
general statements about the structure and degeneracy of the entanglement spectrum 
for blocks of any length. 



2. Model and definitions 



2.1. Quantum algebra 

Let us denote states of a spin- S at site i by \S, m)i. Here m G {— S, — S + 1, . . . , S} 
is the magnetic quantum number denoting the z-component of the spin. The label 
S of the state \S, m)i is invariant under the action of g-deformed angular momentum 
operators satisfying the SU g (2) quantum group algebra [2] 



\2Jf 



[J?,J?] = ±J? 



2/2 



</2 



,1/2 



-1/2 



(3) 



This algebra has two different unitary representations for positive real q and complex 
q on the unit circle [41]. In this paper we will consider the former case where q G R + . 
The resulting algebra is invariant under the transformation q — > g -1 so that we 
consider further q G (0,1]. The usual SU(2) algebra is recovered at the isotropic 
point q = 1, while full deformation occurs in the limit q 
be used extensively below. 

The analogue of total angular momentum, J to t = J i 
of operators through the definition of the coproduct 



0. The g-number [n] will 
J2 is realized at the level 



J tot 
u tot 

TZ,± 



-Jl/2 



j£ + Jf 



q - ^ j 2 
Ii ® J 2 Z + J* ® I 2 



(4) 
(5) 



The operators J^ t satisfy the quantum group algebra (3). A (2 J + l)-dimensional 
irreducible representation of Jtot is therefore spanned by the states 

h h J' 



\J,m)= 



mi 7722 fn 



|jl, 777l) ® \ j2,m 2 ), 



which satisfies 



^tttl^ m ) 

Jtot\J, m ) 

These equations define 



= y/[JTm][J±m+l] \ J,m±l), 

- 771 I J, 77l) . 

the q-deformed Clebsch-Gordan 



(g-CG) 



(6) 



(7) 

(8) 

coefficients 



J K L 
rrij rrik mi 



up to a phase [41-43]. 



The q-CG coefficients are components of a unitary matrix (a change of basis 
matrix) and may be chosen to be real. These coefficients vanish if the triangle relation 
\j\ — 32\ ^ J ^ jl + J2 and selection rule mi + TO2 = to are not satisfied (angular 
momentum conservation). Throughout this paper, the summation indices mi (lower 
row of g-CG symbols) are understood to run over all values compatible with the 
corresponding quantum number ji (upper row of g-CG symbols). For example, in 
(6) we sum over m, G + 1 , . - . , Ji } ■ Some identities involving the g-CG 

coefficients that we use in the following derivations are collected in Appendix A. 
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Figure 1. Diagrammatic representations of matrices used in the MPS description 
of the VBSq(S') state (upper row). The g-deformed Clebsch-Gordan coefficients 
(lower left) are important in the diagonalization of the transfer matrix. The 
matrix p; (lower right) is related to the reduced density matrix p; of I sequential 
spins in a chain of L sites by p; = p; /trG L . 



2.2. Matrix product representation 

Some of the objects we describe here are conveniently represented as diagrams 
(figure 1). The manipulation of these diagrams has been useful in the study of 
entanglement and correlation functions in matrix product states [27,44]. 

Let us now write down the MPS representation of the VBSg(S) state. For a 
periodic chain of L spins we have 

|VBS g (5))=tr(g 1 .g 2 .....g L ), (9) 

where gi are (S + 1) x (S + 1) matrices. The trace here is done over the auxiliary 
matrix space. The elements of gi and its dual g^ arc state vectors: 



(gt)o6 = X! 

m 

(Si)ab = ^2 



S S o 

a —b m 



S S c 
2 2° 

a —b m 



(-l)V b/2 \S,m)i 

i 

(-1)V /2 (S,m\i 



(10) 
(11) 



The |VBS g (S 1 )) state (9) is annihilated by the g-deformed AKLT Hamiltonian (periodic 
boundary conditions) 

L L 2S 

n = J2 h u+i = J2 E n s(M + i), (12) 



i=l s=S+l 



where II s (i,i + 1) is a projector onto the subspace spanned by the g-deformed s- 
multiplct formed by spins at i and i + 1. To prove this, we look at the overlap between 
the two states 



(»-gi+i)ao= E (-ir^ (b - c) 



\J,m) = 



brn'm 

' jl .72 J 

mi nT-2 tti 



r s 

2 


s 

2 


S ' 




- s 

2 


| s- 


a 


b 


to' 


i 


-b 


— C TO 



\ji,mi)i ® \j2,m 2 ) 



i+l- 



\S,m') t ®\S,m) l+1 , (13) 



(14) 
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Figure 2. Diagrams representing the contraction of the transfer matrix G with 
a g-deformed Clebsch-Gordan coefficient. The g-deformed .F-matrix is defined 
according to the upper diagram in which the leg labeled by b' is shifted. The 
lower diagram represents the eigenvalue equation (19). In these diagrams internal 
lowercase indices are summed over. 
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Since the states 777) i} are orthonormal to each other, we obtain 
(J,m\(Si-Si+i)ao= E C" 1 ) 6 " 6 ? 

bmirri2 

Using an identity (B.2) derived in Appendix B yields 
(J,m|(g,.g 4+1 ) ac = (-l) c - s /V c/2 



'[25+1] 



[S+l] 



F q [S^J^; -^S] 



s S j 

2 2 J 
a — c 777 



.(15) 



(16) 



The elements F q [DBJC;NK] of the g-deformed F- matrix are defined diagrammati- 
cally in figure 2. The q-CG coefficient in the overlap (16) vanishes when J > ^ + -| 
proving that |VBS 9 (5)) = 0. Furthermore, has nonnegative eigenvalues 

because it is a sum of projectors. Thus, |VBS 9 (S')) is a ground state of the Hamiltonian 
(12). 



2.3. Transfer matrix 

We now construct a transfer matrix G that is defined in terms of g and g by 
{G)aa'-,bb' = (g)ah (g)a'fc' ■ Explicitly, its elements are 



(G) aa , 



bb' 



E 



r s 

2 


f & 




- s 

2 


s 
2 


S ' 




—b m' 




a' 


-b' 


m' 



(-l) h +V (6+6 ' )/2 . (17) 



This transfer matrix is an important object that appears in the calculation of 
correlation functions and the reduced density matrix of MPS. Let us diagonalizc 
this matrix through an approach based on the g-deformed J 71 - matrix (Appendix B). 
As depicted in figure 2, we construct the eigenvalue equation (G) aa i e^ = Xe aa ' 
using the q-CG coefficients as an ansatz for the elements of the eigenvector e aa i. The 
resulting equation is 



(G) aa '-bb'ebb' 



E 



biz 



r s 
2 


f & 




^ s 

2 


j 


s - 
2 




- 8 

2 


s 
2 


S ' 


a 


—b m' 




_ b 


m 


v. 


Q 


a' 


-b' 


ml 



^_ 1)b+b > q - {b+b >)/2_ (lg) 



Transposing columns in the third q-CG coefficient in order to match the identity (B.2) 
leads to 
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S 1 
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ml 




—6 m 


b'_ 




in' 


V 


a'_ 



\ ' r c I 11 r 1l D 2 2-J' 2 2J 



S „■ Si 
2 J 2 

a m a' 



-1)" 



= Ae„ 



"V /2 i 



[25 + 1] 
[5 + 1] 



(19) 



[5+1] 

We see that the elements of the eigenvectors of G are e aa ' = (ejm)aa' = 
A suitable similarity transformation on G gives the orthonormal set 



S/2 j S/2 
a m a' 



of eigenvectors 



^■jm)aa' 



s s ■ 

2 2-1 



(20) 



To obtain this set we transposed the middle and last rows of (e jm ) aa ' and considered 
the orthogonality relation (A. 2). The eigenvalues associated with these eigenvectors 
are 



A; 



( ir g [ 25 + 1 ] p raSSj.SS] 
y L ) _|_ -n r ll J 2 2J> 2 21 > 



= (-iy +S [2S+l] 



K s 

2 2 
2 2 



(21) 
(22) 



where the g-deformed 6j symbol in the second line is defined in (B.4). The eigenvalue 
Xj is (2j + l)-fold degenerate and < j < 5. The absolute value of Xj decreases with 
increasing j. These expressions match the results of [35, 36] (except for a multiplicative 
constant). 



2. 4- Reduced Density Matrix 

In this subsection we calculate the reduced density matrix pt of I sequential spins in 
a chain of infinite length L — > oo. Using the formalism developed in [5] for matrix 
product states, we obtain 

a'b' 



trG 1 



Integer powers of the transfer matrix G may be written as 



G' 



(24) 



where P, is a projection matrix onto the subspace spanned by the eigenvectors e.j m 
of G. Since the dominant (largest magnitude) eigenvalue of G is Aj = o> large integer 
powers of G simplify to G" — > AqPo as n — > oo. Thus, in the limit of infinite chains 
L — > oo the reduced density matrix (23) simplifies to 

5 X e 

{Pt)ab-cd = ^](Po)aa';6h' ^ (Pj )a'c;b'd- (25) 



ia' ;bb' ^ ^ 
a'b' j=0 ^0 

Constructing the projectors P, from the eigenvectors of G yields 

„-(o+6)/2/_^\o+6+S s 



£ )ab;cd 



[5 + 1] 



■ E 

j—0 m— — j 



E 
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m 



(26) 
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t_l\a+S/2 q -a/2 



r s s o 

2 2 u 

—a a 



(27) 



We can further express the reduced density matrix as a sum of tensor products by 
defining the (5 + 1) x (5 + 1) matrix 

f_]\a+S/2q-a/2 



S S ■ 
2 2 J 



S-n 



(28) 



The Kroncckcr delta here enforces the triangle relation. Finally, making the necessary 
substitutions gives 

s X £ j 

Pe = J^TT Y] Qjm^Qjw (29) 



j—0 ^ m——j 



3. Double scaling limit 

3.1. Eigenvalues of reduced density matrix 

In the double scaling limit, we consider infinitely long blocks and take I —> oo. The 
reduced density matrix pe (29) simplifies to a tensor product of diagonal matrices 
Poo = Qoo ® Qoo- Explicitly, we have: 



ib\cd 



[5 + 1] 



r S S 


0" 




2 2 u 


2 2 




—a a 





g 


-6 6 



(30) 



Remembering that — 5/2 < a,b,c,d< 5/2, with integer steps, we arrive at the final 
expression for the reduced density matrix, 

q -(a+b) 

(Poo)ab;cd = rn , -,„ SacO~bd- (31) 

[5 + iy 

From this expression we can compute all eigenvalues of p x . For example, in the case 
of a g-deformed spin-2 VBS state we have 

q- 1 \ / q- 1 
10® 10). (32) 
q J \ q 



1 



Poo = 



(q + l + q- 1 ) 2 

The nine eigenvalues of this matrix are proportional to {g 2 , g, q, 1, 1, 1, q~ 2 } 



3.2. Entanglement spectrum and entropy 

We now write poo = e~@ Hc / tr e~@ Hc , where H c is an effective Hamiltonian and 
1/(3 an effective temperature. The eigenvalues of the Hamiltonian H c constitute the 
entanglement spectrum of the block [10]. The tensor product form of p^ (31) yields 
the simple paramagnetic model 

-0H e = -P{HP + iff)) = ph{S{ + 5f), S 2 = f(f + 1), (33) 
/3h=\]nq\. (34) 

Here h is the magnitude of an effective magnetic field along the z-axis, while Sj are 
spin-5/2 operators of the undeformed SU(2) algebra. We can thus identify |lng| as the 
ratio h/T e between the magnitude of the magnetic field and effective temperature T e . 
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Figure 3. The Renyi entropy Sr(c?) of a g-deformed spin-S VBS state vanishes 
in the limit q — > + . At the isotropic point q = 1, large blocks are maximally 
entangled (a) = 2 ln(S + 1) . The von Neumann entropy is obtained in the limit 
a — > 1 (bold line). 



(i) 

We observe that the spectrum of He consists of S+ 1 equidistant energy levels. Thus, 
in the limit S — > oo the entanglement spectrum of the block is equal to the energy 
spectrum of two harmonic oscillators with frequency uj (with an S'-dependent energy 
shift). This frequency is related to the deformation parameter through |lng| = cu/T e . 

In this effective picture, the isotropic case q = 1 corresponds to infinite 
temperature or zero field strength. The block is therefore maximally mixed. The 
reduced density matrix has (S + l) 2 nonzero identical eigenvalues. In the opposite 
limit q — > + the effective model corresponds to zero temperature or infinite field 
magnitude. Hence, the block is in a pure state with zero entanglement. 

We use the eigenvalues of the reduced density matrix (31) to compute the Renyi 
entropy 

HS+V) 



Sr(cx) 



<1 




(35) 



This is an exact expression in the double scaling limit (infinite block). Taking the 
limit a — > 1 gives the von Neumann entropy 



{1 _1 S + l S + l ^ 

K +q \ -(S + i) q s 2 +1 +q si ing, 
qi - q 2 q 2 - q 2 J 

{ S + l S + l ~V (I _l S + l s + i 1 

q \- q ~: 1+ {^^-(g+D ^^^ Ug- (36) 
q2-q 2 J [q2-q 2 q — _ q —J 

In order to recover previous results for the spin-1 case [5], we have to rescale q — > q 2 
in (35) and (36). This transformation is necessary because of the different definition 
(3) used here for the deformation parameter q. 

These entanglement entropies are graphed in figure 3 as functions of the parameter 
q for the cases S = 2, 10. At the isotropic point q = 1, the entanglement entropies 
simplify to 

S R (a) = 5 vN = 21n(5 + 1), q = 1. (37) 
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We thus recover previous results [16, 17] for isotropic spin-5 VBS states. In the limit 
of full deformation q — > + the entanglement entropy for any spin 5* vanishes. 

Finally, let us consider the case of very high spin at fixed < q < 1. Taking the 
limit S — !• oo in (35) and (36) gives 

2 

SB,(a) = - In 



S y 



1 - a 
2 ml 



( g -l/2_ g l/2)< 
q -a/2 _ q a/2 



1 



1/2 



-1/2 



lng, 



(38) 



S -> oo. (39) 



q -i/2_ q i/2j ^ g i/2_ g -i/2^ 

We find that the entanglement entropy is bounded for any g-deformed AKLT chain of 
arbitrary spin S. It diverges only at the isotropic point q = 1. 



4. Finite-size corrections 

We now look at the case of finite blocks in an infinite chain. The reduced density 
matrix (29) of a block of £ spins may be written as 



Pi 



An 



3 

E 



Qjm ® Qj n 



Let us express this matrix in the basis of the following vectors: 



{XjAl)ab 

with corresponding dual 



(-l)-(J+b) q b/2 

y/[2J+l] 



S/2 S/2 J 
a -b M 



(vjmU = (-iy j+b q- b/2 vw+n 



S/2 S/2 J 
a -b M 



This choice of basis allows us to use identity (B.2) and get 
{ _ iy >-j 



(Pi) 



J'M'JM 



[5 + 1F 




A 



s si 



± J ^£ r 1 l J 2 2 ° ' 2 2 



S/2 S/2 J 
a -6 M 



(40) 



(41) 



(42) 



(43) 



We find that the reduced density matrix can be decomposed into sectors labelled by 
the quantum number M £ [-S, S]. Each sector is represented by an (S + 1 — \M\) x 
(5+1 — \M |) matrix. These matrices are diagonal in the {vjm} basis only when q = 1. 
We treat this isotropic case analytically in section 4.1. For the general case < q < 1, 
we calculate the eigenvalues of pi using first-order perturbation theory and compare 
this approximation with numerical results in section 4.2. 



4-1- Isotropic case 

When q = 1, using the orthogonality condition (A. 2) in (43) yields the exact (2 J + 1)- 
fold degenerate eigenvalues 



Pjm 



3=1 
\J+S s 



A 



(s + iy 



3 = 1 u 



/f!Ji 



(44) 



(45) 
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with < J < S and M G {—J, — J + 1, . . . , J}. For instance, taking S — 2 gives the 
exact eigenvalues 

Poo = 5 (1 + 3(-2)-' + 5(10)-') , (degeneracy 1), (46) 

pi^ = |(l + f(-2)^-|(10)- / ), (degeneracy 3), (47) 

P2M = H 1 -f(- 2 )"' + 5( 10 )"')' (degeneracy 5). (48) 

The formula (45) reproduces the results of [16, 17] for undeformed spin-S* AKLT 
chains obtained from the Schwinger boson representation of the VBS state. Our 
approach, however, emphasizes the role of 6j symbols in determining finite-size effects 
on entanglement in these states. Additionally, this result solves a recursive formula in 
[16, 17] for the coefficients in the sums for the eigenvalues pjm- 

The leading finite-size correction to the eigenvalue pjm is proportional to the 
exponential factor (Ai/Ao) f = (— l) c c~ £ / ? . Using the formula 

Jsffl (SI) 2 

\iff J, (s-j)\(s+j + i)V m 

gives the characteristic length of decay £ = 1/ In ((5 + 2)/ 5). This length is equal to 
the correlation length of the spin-spin correlation functions in the spin-S VBS state 
[36]. 

Let us construct an effective Hamiltonian for long blocks 1 <C I < oo in the 
isotropic case. Considering only the leading-order correction to the eigenvalues (45) 
gives 

"- " <sTTF j 1 - s(sT2) (sTi)'* 2 ' 71 ' + " - S{S + 2 » } ■ (50) 

Since the reduced density matrix is diagonal in the {vjm} basis, we can write the 
effective Hamiltonian H c as 



-pH B w In {I - 



S(S + 2) 

12 



S(S + 2) 



(^)'(2J(J + l)-S(S + 2)) 
(^2) xi{J(J + l)-5(f + 1)}. 



(51) 



This expression is valid for 35^(S + 2)~ l <C 1. If we define an undeformed spin-S 
operator J = Si + as the sum of two spin--| operators Si and S; on the block 
boundaries, we obtain the Heisenberg model 

^ = 7(^-1)^.8,, 7(^)= (52) 

We can identify T c = l/ r y(S,£) as an effective temperature that depends on 
the length of the block. The double scaling limit £ —> 00 therefore corresponds to 
a maximally mixed state (infinite temperature). In this interpretation, we observe 
that the sign of the coupling strength changes with block length (alternation between 
ferromagnetic and antiferromagnetic interactions). This implies that the dominant 
eigenvalue of the reduced density matrix alternates between the poo singlet (even £) 
and psm multiplet (odd £). 
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Let us now consider the entanglement entropy of a block consisting of a single 
spin (t = 1) for the case q = 1. The eigenvalues of the reduced density matrix may be 
written as 

{ _ iy+ s 



Pju{i = l) = 



1 



+ 



(S+l) 2 ' [SS/2S/2' 
\ S/2 S/2 

Making use of the identity (49) and 
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gives the desired result 




S. 



SJ 



25 + 1' 



PJM(e = i) = 



q=l. 



(53) 

(54) 
(55) 

(56) 



25+1' 

Thus, the single-site reduced density matrix has (25 + 1) nonzero identical eigenvalues. 
This result proves that the block is a uniform mixture of the (25 + 1) states of a 
single spin-5 as expected. The entanglement entropy in this case is Sr. (a) = 5 v n = 
ln(25+l). 

For long blocks satisfying I 3> £, the leading nonvanishing correction to the 
entanglement entropy is proportional to (Ai/Ao) 2f . The approximate Rcnyi entropy 
in this case is 

21 



S R (a)«21n(S + l)- 



3 a 



5 



5(5 + l)(5 + 2). 



(57) 



2 \S + 2 j 

Finite-size corrections to the von Neumann entropy can be obtained from this result 
by taking the limit a — Y 1. 



4-2. Anisotropic case 

For arbitrary values of q, the dominant characteristic length of finite-size corrections 
generalizes to £ = 1/ ln([5 + 2]/[5]). That is, we have |Ai/A | = [5]/ [5 + 2] < 1. For 
blocks of length I ^> £, we may therefore approximate the reduced density matrix as 



pi ~ Qoo ® Qoo + 



t 1 

i E 

m=-l 



Qlm (8 Ql 



(58) 



We have already determined that Q o is diagonal with nondegenerate eigenvalues (28). 
This means that first-order perturbation theory within each sector of the preceding 
equation involves only the diagonal elements of Qi m . From (28) we know that only 
Qio has nonzero diagonal elements and hence we obtain the approximate eigenvalues 



-(M+f) 



Pfj.v — Pull 



[5 + 1] 



ll + [3]^ 



5 1 O 
2 2 

jJL n 



1 1 s 

2 2 

V V 



(59) 



The labels fi and v are quantum numbers that run from — ^ to ^ with integer steps. 
The second term in (59) involving the q-CG coefficients may be evaluated explicitly 
with the identity [41] 



s 1 s 

2 2 
[1 \l 



y/[S][S + 2] 



| 9 |(l+S/2) 



I' 



(60) 
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Figure 4. The eigenvalues of the reduced density matrix of a block of I spins in a 
spin-2 VBS 9 state (solid and dashed blue lines) are compared to the perturbation 
result (red dotted lines). Solid blue lines denote nondegenerate eigenvalues while 
dashed blue lines denote doubly degenerate ones. The dominant eigenvalue 
approaches unity as q — > 0. For 1 = 1 four eigenvalues are zero for all q. 

These approximate eigenvalues are compared to exact numerical results for the spin-2 
case in figure 4. We observe a rapid improvement in the accuracy of the perturbation 
result with increasing block length £. Furthermore, these numerical results reveal how 
g-deformation modifies the degeneracy of the entanglement spectrum by breaking the 
multiplet structure present in the isotropic case. 

5. Conclusions 

We exactly calculated the reduced density matrix of g-deformed VBS states with 
arbitrary integer spin- S in the double scaling limit. We discovered that the 
entanglement spectrum corresponds to a thermal ensemble of two spin-S'/2's in a 
uniform magnetic field h. The deformation parameter q enters this picture as the 
effective field divided by the temperature lnq = h/T c . In this infinite block limit, we 
also derived exact expressions for the Renyi and von Neumann entropies as functions 
of the deformation parameter q and spin S. We used this result to demonstrate that 
entanglement in a VBS (? (S') state is monotonically reduced by increasing anisotropy. 

Furthermore, we obtained the exact reduced density matrix of finite blocks in an 
infinite chain. We diagonalized this matrix for the isotropic case q = 1 and obtained 
its exact spectrum in terms of 6j symbols. We found that degenerate eigenvalues 
of the reduced density matrix are grouped into multiplcts. For long but finite 
blocks, we constructed an approximate effective Hamiltonian in the undeformed case. 
This effective model consists of two spin-5/2's with a Heisenberg interaction. The 
coupling parameter of this Heisenberg Hamiltonian alternates between ferromagnetic 
and antiferromagnetic depending on the parity of £. 

We also studied the general case of arbitrary deformation q and block length £. 
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We made approximations for the eigenvalues of the reduced density matrix of long 
blocks using first-order perturbation theory. Finally, we numerically investigated the 
spectrum of the reduced density matrix of finite blocks in a deformed spin-2 VBSg(5) 
state. In this case we discovered that q-deformation partially breaks the degeneracy 
of eigenvalues within each multiplet. 

We have mentioned in the introduction that generalized VBS states have been 
constructed with various symmetries. The dimension of the MPS representation of 
these states depends on the particular representation chosen for the given symmetry 
group. Let us consider a D x D MPS. On general grounds one can state that the rank 
of the reduced density matrix is less than or equal to D 2 regardless of the number of 
spins in the block [5]. Thus, the entanglement entropy is always bounded according 
to S < 2\nD. For example, the isotropic SU(V) VBS state is described by an N x N 
g matrix. This entropy bound is saturated in the double scaling limit where it is 
known that the entanglement entropy is S — 21nV [16,17]. In this case, we expect 
anisotropy by g-deformation to reduce entanglement entropy in the SU(V) VBS state 
as we proved here (figure 3). 
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Appendix A. Identities for g-CG coefficients 

Among the key properties of the q-CG coefficients that we use above are the 
orthogonality relations 
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' jl 02 


j~ 


m\ ?7i2 m 


i 




m 



jl 32 J 
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3i 32 J' 

mi 1^2 TTl' 



(columns), (A.l) 
= Sjj>S mm ', (rows). (A.2) 



We also make much use of the following identities involving column transpositions: 

= (_l)Jl- J + m 2g- m 2/2 / [2^+1] 



Jl 32 J 

mi mi m 

ji 32 J 
mi m-2 m 

ji 32 J 
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(A.3) 
(A.4) 

(A.5) 



Appendix B. g-deformed F-matrix and 6j symbols 



The equation for the lower diagram given in figure 2 reads 
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Using the identity (A. 3) in the righthand side of (B.l) and applying the orthogonality 
condition (A. 2) to evaluate the sum gives 
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(B.2) 
(B.3) 



Here ^ s the 9-deformed 6j symbol. It is related to the elements of the 

g-dcformcd i^-matrix by [41] 



(DB A\ 


'AC J' 


\CJK] q 
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F q [DBJC; AK] = (-1) D + B + J + C ^\2K + 1}[2A + 1] j £ * £ } 



(B.4) 
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